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ABSTRACT 

The radiation from afterglows of gamma-ray bursts is generated in the collisionless plasma shock interface 
between a relativistic outflow and a quiescent circum-burst medium. The two main ingredients responsible for 
the radiation are high-energy, non-thermal electrons and a strong magnetic field. In this Letter we present, for 
the first time, synthetic spectra extracted directly from first principles particle-in-cell simulations of relativist 
collisionless plasma shocks. The spectra are generated by a numerical Fourier transformation of the electri- 
cal far-field from each of a large number of particles, sampled directly from the particle-in-cell simulations. 
Both the electromagnetic field and the non-thermal particle acceleration are self-consistent products of the 
Weibel two-stream instability. We find that the radiation spectrum from a F = 15 shock simulation show great 
resemblance with observed GRB spectra — we compare specifically with that of GRB 000301C. 
Subject headings: radiation, acceleration of particles, gamma rays: bursts, shock waves, instabilities, magnetic 
fields, plasmas 
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1. INTRODUCTION 

Radiation from gamma-ray burst (GRB) afterglows is gen- 
erally believed to be synchrotron radiation. The main ingre- 
dients in this type of radiation are a strong magnetic field 
and a power-law distribution of high-energy electrons. Multi- 
wavelength observations are more or less well fitted by hid- 
ing our ignorance about the microphysical details of these in- 
gredients in the parameters: eg and describing the amount 
of the total shock energy that is deposited in magnetic field 
and electrons respectively, and p; the slope of the non-thermal 
electron distribution, N('-f)d^ oc 7"^ where 7 = [1 -(v/c)^]"'/^ 
is the relativistic Lorentz factor. However, the true origin and 
nature of these two ingredients remain unexplained. 

Regarding the magnetic field, observations indicate that 
the typical value of th e equipartition parameter in a fterglows 
is ffi = 0.0001-0.1 ("Panaitescu & Kumai"' '2002', 'Waxman' 
[1997; Wijers & Galama 1999; Yost et al. 2003). This is many 
orders of magnitude larger than what can be achieved by 
shock compression of the interstellar medium. Even if the 
field is injected from the progenitor it is hard to maintain 
such a high level of magnetic field in the shock region. 
lOruzinov & Waxman (1999) and Medvedev & Loeb ( 199^ 
suggested that the Weibel two-stream instability can gen- 
erate a strong turbulent magnetic field in the shock region 
iWeibel| | 1959h. This has been confirmed by nu merical simu- 
lations ("Frederiksen et al."2004*, 'Hededaf2005'; 'Hededal et al.' 
pb04; Hededal & Nishikawa 2005; Nishikawa et al. 2003; 
Nishikawa et al. 2005b; Silva et al. 2003). It has furthermore 
been suggested that the radiation from non-thermal particles 
in this turbulent magnetic field c an be described as 'jitter-' or 
'diffuse synchrotron' radiation JFleishmanI l2005t iMedvedevI 
12000,2005). 

Non-thermal acceleration of particles is most commonly 
described as Fermi acceleration (e.g. N iemiec & Ostrowskil 
and references therein). Fermi acceleration relies heav- 
ily on assumptions about the nature of the magnetic field 
in the upstream/downstream region. The Fermi acceleration 
mechanism also faces several observational problems: E.g. 
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in the Crab Nebula, the low energy electrons have a power- 
law distribution spectrum 1.1 < p < 1.3, much lowe r than 
the "universal" p = 2.2-2 A dWeiler & Panagialll978h . Ad- 
ditionally, with Fermi acceleration, one expects the presence 
of an X-ray halo around the shock, but such a halo is not seen 
in Chandra observations of SN 1006 ( Lo ng_et_aL, 2Q03 ) . Re- 
cently Baring & Braby (2004) found that particle distribution 
functions (PDFs) inferred from GRB observations are in con- 
tradistinction with standard acceleration mechanisms such as 
diffusive Fermi acceleration. We stress that the Fermi acceler- 
ation mechanism in principle may work, but that it has never 
been explicitly demonstrated to dominate, or even to be sig- 
nificant, in first principles simulations. 

We are entering an era where very large simulations with 
three-dimensional particle-in-cell (PIC) codes are becoming 
computationally affordable. Such codes work from first prin- 
ciples and can be considered as numerical experiments. 3D 
PIC experiments of relativistic collisionless plasma shocks 
have b een performed by iFrederiksen et al. (2004); HededaJ 
(2005'); "Hededal et al. (2004); Hededal & Nishikawa (2005); 
Nishikawa et al. (2003); Nishikawa et al. (2005b); Silva et al] 
"^003). These simulations have shown that the two-stream- 
instability generated magnetic field is highly turbulent and 
varies greatly through the shock region. In the non-linear 
stage of the instability, filaments merge and can become as 
strong as < 0.1. Recently, Hededal et al. (2004) have fur- 
ther shown that the Weibel two-stream instability also in- 
volves non-thermal acceleration of the electrons. This new 
acceleration mechanism differs from Fermi acceleration in the 
sense that it is local and instantaneous. Electrons are trapped 
in the potential of the generated current filaments and oscil- 
late with strong acceleration / re-acceleration. The numeri- 
cal experiments have also shown that the collisionless shock 
transiti on zone is at least of the order hundreds of ion skin- 
depths (Frederiksen et al.''2004VHededar2005VHededal et aU 
1^04). This is considerably thicker than previously thought 
(Gruzinov 2001). 

In this Letter we present a novel tool that allow us for the 
first time to extract radiation spectra directly from PIC exper- 
iments. In section|2lwe present the radiation tool , in section 
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|3lwe apply the tool to investigate the 3D jitter radiation spec- 
trum, and in section |4] we present synthetic radiation spectra 
from a 3D relativistic shock simulation resembling a GRB af- 
terglow collisionless shock. Finally, in Section [S] we give a 
brief discussion on the results and compare with observations. 

2. SYNTHETIC SPECTRA FROM PIC-CODE SIMULATIONS 

A particle-in-cell (PIC) code simulates plasma on a much 
more fundamental level than MHD. The code works from 
first principles, by solving the Maxwell equations with source 
terms for the electromagnetic fields, together with the rel- 
ativistic equation of motion for a large number of charged 
particles. The electromagnetic fields are discretized on a 
three-dimensional grid whereas the particles are defined with 
continuous positions and momenta within the grid. The 
relativistic PIC code t hat we have used is described in 
iFrederikse n et alJ J2004 and|Hededal et al. (2004). Recently, 
radiative cooling and synthetic spectra generation has been 
added to the code (Hededal 2005). 

As a particle undergoes acceleration it will emit radiation. 
In the PIC code we know this acceleration at each time step 
and hence we can reduce the energy of the particles in accor- 
dance with the Larmor radiation formula. The energy loss is 
highly dependent on the particles Lorentz factor and for rel- 
ativistic particles, most of the radiated energy is beamed into 
a narrow cone with opening angle 1 /7 around the momen- 
tum vector Because of this relativistic beaming, we take an 
approximative approach and subtract the momentum lost to 
radiation along the velocity vector of the particle. We have 
tested this implementation of radiative cooling against the an- 
alytical result for a gyro-orbiting electron in a homogenous 
magnetic field and find excellent agreement (Hededal 2005). 

Given the detailed information of particle positions, veloc- 
ities, and accelerations it is possible to compute the emitted 
radiation spectrum. To do this, we start by adopting the ex- 
pression from Jackson ( 1999) for the retarded electric field 
from a charged particle moving with instantaneous velocity 
/3 = v/c under acceleration /3 = v/c. 
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Here, q is the unit charge, c is the speed of light, eo is the 
vacuum electric permittivity, and R is the distance to an ob- 
server, n is a unit vector that points from the particles retarded 
position towards the observer. In Eq. [0 we have neglected 
a "velocity term", which is justified when the typical length 
scale of the emitter is negligible compared to the distance to 
the observer (Jackson 1999). The radiation spectrum from an 
accelerated charge is then found by Fourier transforming the 
Poynting flux based on Eq.[2in the retarded time frame. This 
gives us the energy dW radiated per unit solid angle dft and 
unit frequency duj 



d^W 



dVlduj 167r3 



nx [(n-/3)x/3] ,.^(,/_n. 
(l-n./3)2 ' 



(2) 

This expression can only be analytically integrated under 
constrains on the morphology of the electromagnetic field. 
For example, in the case of a relativistic electron moving 
in a homogenous magnetic field with zero electric field, 
the integral reduces to the well known synchrotron case 
dRvbicki & Lightman|[T979l) . In the general electromagnetic 



case, no analytical solution exists. However, from the PIC 
experiments we readily have positions, velocities and accel- 
erations for the particles that define the plasma. With this 
information, we can numerically integrate Eq. |2] to obtain a 
synthetic radiation spectrum from each of the particles in a 
plasma. Finally we can add the spectrum of each particle into 
a total spectrum. Adding the spectrum from each particle lin- 
early is valid as long as the phase of each contribution is com- 
pletely uncorrected to the others. 

We have implemented Eq. |2] and tested several scenarios 
where an analytical solution exist. This includes synchrotron 
radiation, bremsstrahlung, and wiggler/undulator radiation. 
In all tests we find that our radiation tool can effectively re- 
produce the theoretical expected solutions. FigureQlshows the 
synthetic spectrum (solid blue) of a single relativistic electron 
gyrating in a homogenous magnetic field. The envelope curve 
(dashed red) shows the analytical solution, which shows good 
agreement with the simulations. The inset in Fig. \l\ shows 
that the peaks in the spectrum fall on the relativistic gyro- 
frequency and its higher harmonics as expected. For more 
details and tests, see jiededa l (2005). 
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Fig. 1 . — The spectram obtained from a single particle gyrating in a mag- 
netic field from the code (solid blue). The dashed red line shows the theo- 
retical synchrotron spectrum. The dashed line indicates the power-law slope 
2/3. The inset in the upper left comer shows that the spectrum consists of a 
discrete set of spikes that are integer overtones of the gyrofrequency, in this 
case LOB = 0.094 (arbitrary units) (dotted vertical lines). The particle has a 
Lorentz factor 7 = 5. All numbers are in simulations units. 



3. EXAMPLE: 3D JITTER RADIATION 

The PIC plasma simulation code, in combination with the 
radiation generation module, provides us with a powerful tool 
for testing various non-linear problems. In this section we 
examine the so-called jitter- or diffuse synchrotron radiation 
spectrum (Fl eishman 2005: MedvedeY.20Q0...2QQ5 ). This type 
of radiation is generated by relativistic particles that are de- 
flected in a turbulent magnetic field, where the deflection an- 
gle is comparable to the relativistic beaming angle. 

We setup PIC experiments with a periodic random magnetic 
field with a power spectrum that follows a power-law distri- 
bution in the Fourier domain PB{k) oc fc^ (Fig.|2j. In this field 
we use the PIC code to trace an isotropic momentum distribu- 
tions of electrons with a Lorentz factor 7 = 3. The strength of 
the magnetic field is set such that the maximum deflection of 
the electrons is smaller than the relativistic beaming angle for 
all Fourier nodes. 
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Fig. 2. — Left column: The Fourier structure of the turbulent magnetic field 
used in the jitter radiation experiments. Here we show three different values 
of the spectral power law index /j: = — 2 (red noise, top), fi = (white noise, 
middle) and fi = 2 (blue noise, bottom). Right column: CoiTesponding spatial 
2D-slices of the magnetic field. 



Fig. 13 shows the resuhing radiation spectrum for different 
values of the slope ji of the magnetic Fourier decomposition 
Psik) oc k^^. For < 0, the high-energy part of the spectrum 
follows a power-law with a slope a ~ /i - 1 independent of 
the electron energy. For /i > 0, the flat spectrum continues to 
higher frequencies than for /i < and the spectrum has a hard 
cut-off at high frequencies. In all cases, the low energy part is 
flat and quite similar to the case of relativistic bremsstrahlung. 

We observe that as the strength of the magnetic field is in- 
creased, so that the typical deflection angle of the electrons 
becomes larger than the relativistic beaming angle, the re- 
sulting spectrum converges through the wiggler domain to the 
synchrotron case Attwood (2000) . 

4. RADIATION FROM COLLISIONLESS SHOCK 

Our goal is to obtain spectra from the PIC shock experi- 
ments that may be compared directly with observations. Even 
though it is beyond the scope of this Letter, this will eventu- 
ally allow us to put constraints on the conditions and physics 
of GRB afterglows such as GRB environments and jet struc- 
ture. 

We use the nu merical PIC experiment studied in 
iHededal et alJ i2004h . Two colliding plasma populations are 
studied in the rest frame of one of the populations (down- 
stream, e.g. representing the shocked ISM medium at the 
head of a GRB jet). In this frame of reference, a less dense 
population (upstream, representing the unshocked interstellar 
medium) is continuously injected at the left hand side bound- 
ary, with a relativistic velocity corresponding to a Lorentz 
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Fig . 3 . — The jitter spectrum from a mono-energetic ensemble of electrons 
(7 = 3) moving in a turbulent magnetic field. The different graphs are for 
simulations with different values of 



factor F = 15. The two populations initially differ in den- 
sity by a factor of 3. We use a computational box with 
125 X 125 X 2000 grid points and a total of 8 x 10** parti- 
cles. The ion rest-frame plasma frequency in the downstream 
medium is ojpi = 0.075, rendering the box 150 ion skin depths 
long. The electron rest-frame plasma frequency is ujpe = 0.3 in 
order to resolve also the microphysics of the electrons. Hence, 
the ion-to-electron mass ratio is nii/me = 16. 

As the two plasma populations interpenetrate, we observe 
how the Weibel two-stream instability is excited and current 
filaments are formed (Frederiksen et al. 2004; Hededal et ap 
12004* 'Medvedev & Loeb 1999; Nishikawa et al. 20031 
Nishikawaet al. 2005b; Silvaet al. 2003; Weibel 1959). The 
current filaments induce a highly intermittent electromagnetic 
field with a strength of up to 10% of equipartition. The nature 
of the magnetic field is turbulent, with a po wer spectrum that 
follows a power-law in the Fourier domain aFrederiksen et alJ 
2004). We note however, that unlike the random fields in 
section |3l there exists phase-correlations in the generated 
magnetic field that give the field a more coherent structure. 
The instability also directly leads to acc eleration of electrons 
through a newly discovered mechanism ('Hededal et al."2004t 
Nishikawa et al. 2005a). Electrons caught within the Debye 
cylinder surrounding the ion current filaments are strongly 
and repeatedly accelerated and decelerated. The acceleration 
is instantaneous and local, and thus differs substantially 
from recursive acceleration mechanisms such as Fermi 
acceleration. The resulting particle distribution function has 
a non-thermal high-energy component, correlated with the 
distribution of the generated electromagnetic filaments. We 
note in passing that a fraction of the ions are scattered back 
into the interstellar medium. Thus, ion Fermi acceleration 
may possibly take place in coUisionless shocks. 

To obtain the total spectrum from the accelerated particles 
in the Weibel-generated field, we have traced 20,000 parti- 
cles from the PIC simulation described above. Placing the 
observer somewhere along the jet direction we have numeri- 
cally integrated Eq.|2for each of these particles. 

Fig. 13 shows the radiated spectrum from the coUisionless 
shock simulations described above. Here, we have rescaled 
the frequency into real space units and taken the relativistic 
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Fig. 4. — The spectrum from a F = 15 plasma shock propagating into an 
ISM-like medium. The spectrum peaks in the far infrared. Below the peak, 
we find a power-law segment with slope P(ui) ex. up-l^ . For frequencies above 
the peak frequency, the spectrum continues into the near infrared/optical 
band, following a power-law with slope P(ijj) oc o;"'', where j3 = 0.7. 



Doppler-shift into account. We find that the spectrum peaks 
in the far infrared. Below the peak, we find a power-law seg- 
ment with slope P(uj) cx lu^^^. This is interesting because it 
is steeper than the limiting 1/3 slope for synchrotron radia- 
tion, and may reconcile theory and observations on the "line 
of death" issue (Preece et al. 1998) of observed GRB spec- 
tra. For frequencies above the peak frequency, the spectrum 
continues into the near infrared / optical band, following a 
power-law with slope P(uj) oc uj~^ with /3 ~ 0.7. 

From the power-law slope in the spectrum with (/3 = 0.7), 
the standard procedure commonly used in the GRB commu- 
nity would indicate that the electrons have a power-law distri- 
bution, with slope p = 2A (determined by solving (p-l)/2 = 
(3). The standard conclusion would be that this is consistent 
with Fermi acceleration. However, we strongly emphasize 
that the electrons that have produced the spectrum in Fig. 
are not accelerated by a recursive mechanism such as Fermi 
acceleration. The electrons are instantaneously accelerated 
and decelerated in the highly complicated electric and mag- 
netic field near the ion current channels under the emission of 
strong radiation. This differs substantially from the iterative 
acceleration in Fermi acceleration. Furthermore, the paths of 
the accelerated shock electrons are closer to a random walk 
than circular, which means that synchrotron radiation is not 
an adequate formalism in the current context. 

5. DISCUSSION AND CONCLUSIONS 

The main source of information we have about gamma-ray 
bursts (GRBs) is from multi-wavelength observations of GRB 
afterglows. The radiation from GRB afterglows is generated 
in collisionless shocks between the external plasma (ISM) and 
the GRB jet. This emphasizes the crucial importance of a full 
understanding of the physics of and the radiation mechanism 
in collisionless shocks. In the current standard model, the ra- 
diation is assumed to be emitted as synchrotron radiation from 
power-law distributed electrons in a magnetic field. Follow- 
ing our lack of knowledge about the plasma-physical details, 
it is common practice to parameterize this ignorance with the 
dimensionless parameters e^, and p {eg and e^, express the 



fraction of the total internal energy that is deposited in mag- 
netic field and electrons, respectively, and p is the slope of the 
supposedly power-law electron momentum distribution func- 
tion). 

Since collisionless shocks are extremely non-linear and 
self-interacting systems, the only way to study their micro- 
physics with any credence is through self-consistent, three- 
dimensional relativistic particle simulations. Using this ap- 
proach we have previously explained the origi n and nature of 
a strong electromagnetic field in the shock ( Frederikse n et alJ 
2004) and been able to identify a new non-thermal electron 
acceleration mechanism (Hededal et al. 2004). These two 
ingredients are direct and unavoidable consequences of the 
Weibel two-stream instability. 

In this Letter we have developed a novel tool that allows us 
to extract radiation spectra directly from the particle-in-cell 
experiments. Rather than applying the synchrotron approx- 
imation (assuming homogenous magnetic field and a single 
electron power-law distributed population) we trace a large 
number of electrons in the experiments, and calculate the ex- 
act radiation emitted from each electrons. The tool has been 
thoroughly tested and successfully reproduces spectra from 
synchrotron radiation, bremsstrahlung and undulator/wiggler 
radiation from small-angle deflections (Hededal 2005). 

The tool has been utilized for a parametric study of 3D jit- 
ter radiation, where a population of electrons radiate as they 
move in a turbulent magnetic field whose power spectrum fol- 
lows a power-law in the Fourier domain Psik) oc A;^. We have 
examined how the slope /i influence on the resulting radia- 
tion spectra. For /i < 0, the high-energy part of the spectrum 
follows a power-law with a slope a ~ /i - 1 independent of 
the electron energy. For /i > 0, the flat spectrum continues to 
higher frequencies than for /i < and the spectrum has a hard 
cut-off for high frequencies. 

We have furthermore examined the radiation from a typi- 
cal GRB afterglow collisionless shock that propagates with 
r = 15 through the interstellar medium. We find that the re- 
sulting radiation spectrum peaks just below lO'^Hz. Above 
this frequency, the spectrum follows a power-law F oc , 
with (3 = 0.7. Below the peak frequency, the spectrum fol- 
lows a power law F oc ly" with a ~ 2/3. This is steeper than 
the standard synchrotron value of 1 /3 and thus more compat- 
ible with observations. Both the slope and the peak is con- 
sistent with observations (e.g. Panaitescu (2001) who finds 
13 = 0.67 ± 0.04, a peak at 3 x 10"Hz and T > 10 for the af- 
terglow of GRB 000301C after five days). 

Finally, we would like to stress the following interesting 
point: If one hides all the real physical details of the mag- 
netic field in the dimensionless parameters eB, £e, and p the 
conclusion from standard synchrotron radiation theory would 
be that the slope of the electron distribution A^(7) oc 7"'' is 
found by solving -/? = -(;?- l)/2 = -0.70 79 = 2.4 (which 
is consistent with standard analysis of the GRB 00030 IC af- 
terglow; Panaitescu 2001). This would appear to be con- 
sistent with Fermi acceleration and the supposedly universal 
p ~ 2.2 ±0.2. But a closer look reveals that the accelera- 
tion mechanism is not Fermi acceleration but rather a natural 
consequence of the Weibel two-stream instability. The elec- 
trons are instantaneously accelerated and decelerated in the 
highly complicated electric and magnetic field near the ion 
current channels. Strong radiation is produced in this process. 
The electron distribution behind the radiation is a mixture of 
a thermal component for low energies and a power-law com - 
ponent with p = 2.1 for high energies jHededal et ailll?)?)! . 
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The paths of the electrons are more random than circular and 
the electron distribution function varies with shock depth, as 
does the magnetic topology and strength. 

The outcome of the exercise is thus that it is indeed pos- 
sible, from first principles, to produce synthetic spectra very 
similar to the ones that are observed in GRBs, and that the fit 
of such spectra with the standard model assumptions may be 
more or less fortuitous. 
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